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We study linear electron transport through a molecular wire sandwiched be- 
tween nanotube leads. We show that the presence of such electrodes strongly 
influences the calculated conductance. We find that depending on the quality and 
geometry of the contacts between the molecule and the tubular reservoirs, linear 
transport can be tuned between an effective Newns spectral behavior and a more 
structured one. The latter strongly depends on the topology of the leads. We 
also provide analytical evidence for an anomalous behavior of the conductance as 
a function of the contact strength. 

PACS numbers: 73.50.-h, 73.61.Wp, 85.65. +h 
The success of semiconductor industry is evidently represented by the em- 
pirical exponential "law" for the transistors density as a function of time 
(known as Moore's law [1]). Its extrapolation would predict an atomic size 
gate length in less than two decades. The length domain in between is 
manifestly quantum mechanically dominated: here conventional methods 
adopted to characterize electronic devices [2] have to be updated to more 
sophisticated ones, including a microscopic treatment of electronics at the 
molecular scale [3]. 

The interest in the basic mechanisms of conduction across molecular junc- 
tions bridging metallic pads, already object of scanning tunneling micro- 
scope oriented research, has been intensified by recent experimental achieve- 
ments. For example, I—V characteristics of a benzene-l,4-dithiol ring [4], 
and a poly(G)-poly(C) 30 base pair long double stranded DNA [5] have 
been reported. In a parallel development the use of networks of carbon 
nanotubes (CNT) has been the focus of intense experimental and theoreti- 
cal activity as another promising direction for building blocks of molecular 
circuits [6]. CNT are known to exhibit a wealth of novel properties depend- 
ing on their nano— metrical diameter, orientation of graphene roll up, which 
is conventionally characterized by means of the chiral couple (n, m), and 
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whether they consist of a single (single-wall) or many (multi-wall) cylindri- 
cal surfaces [7]. Recent experiments also confirm the possibility of producing 
heterojunctions made of nanotubes [8]. When such junctions involve differ- 
ent materials [9], the characterization of contacts becomes a fundamental 
issue. This could be the case when CNT are attached to a single molecule 
or a molecular cluster with a privileged direction along the current flow, 
namely a molecular wire. CNT have also been recently employed in similar 
configurations for enhancing the resolution of scanning probe tips [10]. In 
such systems the traditional picture of electron transfer between continuum 
state donor and acceptor species has to be reconsidered in view of the fact 
that the electrodes are mesoscopic and the molecule can bear an electric 
current [11]. 

In this paper we drop the conventional assumption of a continuum of (free 
or quasi-free) lead states. This is a reasonable assumption for mimicking 
the presence of large reservoirs provided by bulky electrodes, but it may be 
inadequate when the lead lateral dimensions are of the order of the bridged 
molecule. This is indeed the case for CNT. 

In what follows, we assume a homogeneous tight binding chain as a model 
molecular wire to isolate contact effects. Electrodes are taken as either 
square lattice tubes (SLT), obtained by imposing periodic boundary con- 
ditions on the longitudinal cuts parallel to the lattice bonds, or CNT. The 
latter consist of a rolled stripe of a graphene honeycomb lattice, also paral- 
lel to the lattice bonds. This configuration generates armchair single-wall 
CNT, (1,1). The consideration of SLT delivers additional analytical insight 
to the numerically studied model of CNT [12]. 

The electronic Hamiltonian H = -fftubes + H w { Te + H coup y mg of the system 
reads 

H = E E (% S n a ,< ~ T^nj,)) K> «| (1) 
o=L,R,wire n a ,n' a 

- ^ T |m L )(n wirc = l| - r |m R )(n wirc = iV| + H.c. 

m L <M c m a <M c 

Here, ^ L ' R , ~f mie ) and V are nearest neighbour hopping terms between atoms 
of the left (L) or right (R) leads, molecular bridge, and the bridge/lead 
interface, respectively; e wiro is the on-site or orbital energy of each of the 
re wirc = 1) ■ ■ • j N chain-atoms relative to that of the leads, e L,R . Note that 
nL,R is a two-dimensional coordinate spanning the tube lattice. Summations 
over raL and mR run over interfacial end-atoms of the leads. In general, 
there are M such atomic positions depending on the perimeter of the tubes, 
whereas, the number of hybridization contacts ranges from a single contact 
(SC) M c = 1, to multiple contacts (MC) M c = M. 



acta printed on February 1, 2008 



3 



In order to derive transport properties, we make use of the Landauer the- 
ory [13] which relates the conductance of the system to an independent- 
electron scattering problem [14]. The electron wavefunction is assumed to 
extend coherently across the device and the two-terminal, linear-response 
conductance g at zero temperature is simply proportional to the total trans- 
mittance for injected electrons at the Fermi energy E-p- That is, g = 
(2e 2 /h)T(Ep), where T = Ylj h j R I ^'ljr 1 2 * s a sum over scattering matrix el- 
ements labelled by open channels in the leads jlr- The factor two accounts 
for spin degeneracy. The transmission function can be calculated from the 
knowledge of the molecular energy levels, the nature, and the geometry of 
the contacts. It can be obtained by solving the Lippman-Schwinger equa- 
tion for H and can be written in the following convenient form [11, 15], 

T(E)=AA\(E)A N (E)\G 1N (E)\ 2 , (2) 

where Gin is the Green function element connecting the two ends of the 
TV-atom-molecule and A L ^ is the left (right) lead spectral density. The 
latter is related to the semi-infinite lead Green function matrix Gi ea d and 
is minus the imaginary part of the self-energy 

^nZ™ R = ^2 ( n ^irc\H\n a ) Gi ea(1 (n a ,n' a ) (n' a \H\n wirc ) . (3) 

n a ,n' a 

The calculation of the spectral function Sn^ rc simplifies due to the form of 
the interfacial coupling in our model. Assuming identical electrodes, lead- 
indices are dropped. Hence, the self-energy becomes 

S = lf 12 G lcad (m,m'), 

C m,m'<M c 

where only surface terms enter in the sum over the states in the leads and 
the effective coupling is defined as r e g- = TyJM c . 

In Fig. 1 the spectral function A = — ImS is plotted for both SLT and CNT 
leads in the SC and MC multiple contact configurations. As a function of 
the number of contacts M c the system interpolates two different scenarios. 
In the MC case the spectral density is effectively the spectral density of one- 
dimensional leads, obtained by Newns in his theory of chemisorption [16]. 
Only the channel without modulation in the transverse profile of its wave- 
function contributes to transport [17]. The two-dimensional character of 
the leads enters as an energy shift of 2^ L,R for SLT and of j L,R for CNT, 
yielding an asymmetric density profile with respect to the atom on-site en- 
ergy e L,R . In contrast, the SC spectrum is symmetric and richer due to 
the contribution of all available channels. Additional features characterize 
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Fig. 1. Spectral density plots for square lattice tubes with circumference of M = 
20 atoms (left panels) and (10,10) armchair carbon nanotubes (right). The two 
different scenarios of a single contact, M c = 1 (top), and of multiple contacts, 
M c = M = 20 (bottom), are shown. 



CNT leads. The minima in A (upper right panel of Fig. 1) are responsible 
for antiresonances observed in the conductance spectrum [12]. 
In eq. (2) Gin can be expressed in the following form [18] 

1 1 SI ( £ \ 2 1 



7 wirc GiAr £(N) 7 wirc C(A r -l) V7 wirc / C(N-2) 

where £(iV)/7 wire is the matrix element Guy for the isolated iV-atom molecule. 
It is given in closed form as a function of the normalized energy £ = 
(E - e wirc )/(2 7 wirc ) and N by 

&(A0 = 



N+l ( l N W+l 

1 



£ + V£^l) - [£ - V£ 



In the limit of weak contact coupling the behavior of the Gin element is 
dominated by C(-W) leading to N transmission resonances in the conductance 
of unit height. Nevertheless, if the effective coupling between the molecule 
and the lead is much larger than 7 wire , £(jV — 2) will become the dominant 
term. As a consequence the conductance spectrum is effectively that of an 
(N - 2)-atomic wire [12]. 
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In conclusion, we have pointed out that the conductance of molecular wires 
attached to mesoscopic leads strongly depends on the geometry and dimen- 
sionality of the contacts. In general, a detailed account of the lead states 
must be provided for describing transport. However, in the multiple contact 
mesoscopic limit conductance becomes independent of the topology of the 
tubular electrodes and transport is effectively one-dimensional. In addition, 
we have explicitly demonstrated the previously observed dependence of the 
effective wire length on the coupling strength. 
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